Modeling Effect of ^ r £“ d ffl°tuned 
and Friction ^ p onse 0 f Turbopumps 



Final Report 

Contract No. NAS8-38348 

George C. Marshall s P ac ®^''^ Ministration 
National Aeronaut.es and s P ace * 

ZrsMI space 358,2 


j. h. Griffin and M-T. Yang 

tePartmen '^SeM"“5 

Pittsburgh, PA 1521 


modeling the 

(N cc?r C nF SHROUD CONTACT AND 

EFFECT OF *1 rup MISTUNEO 

PSTCTION DAMPFRS ON TH -_ , 
*«P0NSE 0F T««w;u B PS .na _ ^ 

?SS?nU..-H«noi 79 P 


N95-20795 


Unci as 


January 


G3/37 0091091 


INTRODUCTION AND BACKGROUND 


Originally the goal of this contract was to deve.op a method for analyzing the s romura, 
response of mistuned bladed disk with nonlinear constraints such as shroud contact a, 
friction damping The approach was to used harmonic balance methods of approx, mati g 
nonlinear response, and efficient, receptance methods for represenring 
„ was postulated that the receptance approach could be used to reduce the "umber °f deg e 
of freedom to just those directly associated with the nonlinear constants The 
balance method of deal, ng with the nonlinearities would replace a system of non .near ord.nary 
differential equations with a system of nonlinear algebra, c equauons. Thus, using b 
a Caches the calculation of the response of a mistuned bladed disk, such as the NAS A 
tlpump would be reduced to dta, of solving a relabvely small system of nonlmear algebra, c 

equations 

Several primary objectives had to be achieved in order to meet these objectives. They 
were: 

1, The response of the system with nonlinearities had to be nearly harmonic and have 
unique solutions 

2 The receptance method had to be developed to a point where it accurately predicted 
system response, and, yet was computationally efficient. 

initial studies were done during the first two years of tins contract Jdresring Urn 
,ssue of representing nonlinearities Two types of nonlineant.es were studied the .rnpact 
nXUlriTg shrouds, and friction dampers. It was found tha, for the case of ne^hbonng 
shrouds hitting each other that the problem could have multiple solut.ons, even in the case 
one blade. This meant that it was not feasible to include this type of phys.cs (shroud l.ft- 

or shroud hitting) when trying to simulate the response of bladed d,sks ! ,^-f p "“blem 

of possible solutions that could occur with mult, pie blades would make this proW* 
extremely difficult to solve. It was found in the initial study, that fncuon c°ns.ra,n, on^ 
" a, hand always resulted in unique solutions that cmrld be represented very we 1 using 
hlonl balance methods. Consequently, i, was decided to go ahead with the second par, of 
the project, but include only fric, ion-type nonlineariues ,n our formulation 

A result of this original study was tha, i, was shown that some very jusual dynamic 
response dta, had been observed in a NASA turbopump was very poss,bly due 

contact. This work is documented in [1], 
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A formulation was then developed for reducing the number of degrees of freedom in a 
mistuning analysis Initial work had been done in developing this approach when NASA 
decided that it had to significantly reduce funding for this project. As a consequence, it was 
agreed that the statement of work for this contract would be revised and the deliverables 
changed to the following: 

• Complete the development of a computer code for performing linear mistuning 
analyses with reduced order models developed from a SPAR/EAL finite element model 
of the blade and an annular disk. Reduced order models are models that can be 
explicitly derived from the finite element model and that have significantly fewer 
degrees of freedom than a comparable finite element model, e.g., 10,000 degrees of 
freedom per blade reduced to six. 

• The code will be based on a Monte Carlo type analysis in which the frequencies of the 
blades vary randomly and the forced response is calculated over a range of excitation 
frequencies 

• Document linear mistuning code 

• Prepare a paper on the mistuning analysis of a NASA turbopump for the Advanced 
Earth -to-Orb it Propulsion Technology — 1994 Conference. 

Consequently, the final effort on this contact has focused on developing a linear mistuning 
code and that development and documentation is the focus of this report. A paper on the 
mistuning analysis of a NASA turbopump is included in volume II of the Proceedings of the 
Advanced Earth-to-Orbit Propulsion Technology — 1994 Conference and is cited as 
reference 2. 


DEVELOPMENT AND DOCUMENTATION OF LINEAR MISTUNING 
CODE 

The development and documentation of the linear mistuning code is documented in the 
three appendices of this report They are: 

• Appendix A: Summary Report. This summary report provides a brief summary of 
the theory and application of the linear mistuning code. 

• Appendix B: User Manual. This manual specifies the input requirements for running 
the linear mistuning code 

• Appendix C: Theoretical Manual This manual provides a detailed explanation of the 
theory and application of the linear mistuning code (in effect, this is an expanded 
version of Appendix A). 
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DELIVERY OF LINEAR MISTUNING CODE 

The linear mistuning computer code (LMCC) and the input files used in analyzing the 
NASA turbopump have been transferred to NASA Marshall Silicon Graphics Computer on 
which the researchers where given an account. The example case, the NASA turbopump, ran 
correctly on the NASA computer There was a problem with compatibility in transferring to 
the Silicon Graphics computer at NASA in that the random number generating subroutine did 
not work correctly on the Silicon Graphics if certain optimization options were used during 
its compilation. This problem is documented in the User Manual. 

Thus, LMCC and example files have been successfully delivered to NASA Marshall at 
this time. 


Interesting results from running LMCC are discussed at some length in Appendix A 
Briefly, however there are two key results: 

1 It works well, i.e., it is computationally efficient and accurate. 

2 When the disk is stiff, it predicts very different mistuned response from the mass 
spring models that were used in the past. LMCC shows that bladed disk systems 
can have a large amount of amplitude scatter from blade to blade even when the 
disk is very stiff The old mass/spring models predict that, in this case, the blades 
would act as uncoupled, independent systems (the disk is almost a rigid 
foundation) and that there would be very little amplitude scatter This could be a 
very important result if NASA decides to go to bladed disk systems with stiff 
disks e g if they decide to use integral bladed disks (Experimentalists at the 
engine companies have also observed large blade to blade scatter and reported that 
the mass spring models did not predict this behavior.) 



Several things could be done to make the linear mistuning code more efficient, more 
accurate, or more convenient. The primary ones are: 

1 Modify the input data so that the disk modes could be calculated using cyclic 
symmetric boundai^ conditions in EAL. Currently, LMCC requires that the 
modes of the disk be calculated using a model of the full disk 
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, i-i a Purrentlv only blade and disk 
2 Add aerodynamic coupling between the blades^ C ^ tion was not judge to 
2 md damping are » very shon sttff 

be very important since most NASA 

3 element c^er^Cuneml^i^only^ 0 ^ 1 ^^^^ 50 ^^^'^^^ 

jstjzzzz r « 

man,,er m . asses si„ g various strategies for actively 

Lastly, LMCC could be a „ that its structural ridelity would allow 

controlling bladed disk response The teas re sp„„se of putting actuators a, 

the researcher to quantitattve assess the effect 
various locations on the disk and blades 
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Appendix A: Summary Report 


Executive Summary 


A reduced order approach is introduced that can be used to predict the steady-state 
response of mistuned bladed disks This approach directly takes results from a finite 
element analysis of a tuned system and, based on the assumption of rigid blade base 
motion, constructs a computationally efficient mistuned model with a reduced number of 
degrees of freedom. Based on a comparison of results predicted by different approaches 
it is concluded that: the reduced order model displays structural fidelity comparable to 
that of a finite element model of the entire bladed disk system with significantly improved 
computational efficiency, and under certain circumstances both the finite element model 
and the reduced order model predict quite different response from simple spring-mass 
models. 

Introduction 


The resonant amplitudes of turbine blades tend to be sensitive to minor variations in 
the blades’ properties It is realized that, because of the rotational periodicity of its 
geometry, a bladed disk usually has natural frequencies that are clustered in narrow 
ranges. When the natural frequencies of a system are close together, slight variations in 
the system’s structural properties can cause large changes in its modes, and, 
consequently, its dynamic response. The sensitivity of a bladed disk’s dynamic response 
to small variations in the frequencies of the blades is referred to in the literature as the 
blade mistiming problem and has been studied extensively, for examples refer to Dye and 
Henry (1969), Ewins (1988), Fubunmi (1980), Griffin and Hoosac (1984), or Ottarsson 
and Pierre (1993). It is important to understand mistuning since it can result in large blade 
to blade variations in the vibratory response and the high response blades can fail from 
high cycle fatigue 

Much of the work that has been done in mistuning utilizes spring-mass models to 
represent bladed disks in order to reduce the number of degrees of freedom and to make 
the problem computationally tractable, for examples refer to the previously cited papers 
The model’s parameters, such as the mass and the spring constants, are chosen in an ad 
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hoc manner and one must question the ability of such simple models to accurately 
represent such complex systems While some attempt has been made to corroborate the 
accuracy of spring-mass models by comparing predictions with specific test data, for 
example Griffin (1988), such work is relatively scarce 

Efforts have been made to develop more structurally accurate models for bladed disks 
by using plate elements to represent the disk and beam elements to represent the blades, 
for examples refer to Kaza and Kielb (1984), and Rzadkowski’s two papers in (1994) 
While there can be blade configurations for which the beam representation may be 
adequate, plate, thick shell, and even solid elements are often needed to represent modem 
low aspect ratio blades. The finite element method could be a possible choice to 
accurately model a whole bladed disk, but it is recognized that the time cost and t e 
storage space required to run these programs would be prohibitively high For example 
one could imagine using the Monte-Carlo approach of Griffin and Hoosac (1984) with 
detailed finite element models of the entire mistuned bladed disk. Such an approac 
would involve the analyses of hundreds of mistuned disks in order to determine t 
statistical variations in the blades’ vibratory response. Clearly, such computations are 
currently beyond the capabilities of even super computers and would be hardly suita e 
for use as a design tool Furthermore, because of the extremely large number of degrees o 
freedom involved and the closeness of the natural frequencies, one must question if sue 
results would even be numerically accurate 

The limitations associated with spring-mass and beam models and the direct finite 
element approach motivate us to consider the possibility of developing a new model for 
analyzing mistuned bladed disks. Our goal is to develop a methodology that will direct y 
take the results from a finite element analysis of a tuned system and construct a 
computationally efficient mistuned model with a reduced number of degrees of freedom 
The intent is that the approach will display structural fidelity comparable to a finite 
element model and computational efficiency more comparable to that of a spnng-mass 

model. 



In the study of the steady-state response of complex structures, one widely used 
analytical approach is the receplance method. Bishop and Johnson (1960) The 
receptance method is based on the observation that the dynamic response of every 
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substructure is determined by how it interacts with its environment at its boundaries. If 
the substructure interacts with its environment only at limited areas, it is convenient to 
express the degrees of freedom of the entire substructure in terms of the degrees of 
freedom of its interfaces. The benefit of this method is that when severaJ substructures 
interact with each other, it is only necessary to solve for the degrees of freedom 
associated with the interfaces Once the degrees of freedom of the interfaces are 
determined the response of all substructures and, consequently, the whole structure may 
be calculated To apply the receptance method to the mistuned bladed disk, it is divided 
into two substructures — the disk (Figure 1) and the blades (Figure 2). Modal analysis is 
then used to determine the substructures’ behavior in terms of the degrees of freedom of 
the interfaces. However, the receptance method has two shortcomings: 




1 . The substructures’ modes have to be free at the disk-blade interfaces in order to be 
admissible This works reasonably well for the disk since the blades provide 
relatively little constraint at its rim and, consequently, only a few families of disk 
modes are required to represent its response. However, it is undesirable to use 
the free-free blade model. Because blades generally vibrate close to the clamped- 
free condition, a large number of the free-free modes are needed to achieve a good 
representation of its vibratory response. 


2 . A direct application of the receptance method results in a formulation with a 
reduced number of degrees of freedom. However, depending on the number of 
nodes at the blade interface, the number can still be quite large, especially if solid 
elements are used to model the blade’s neck. Since it is undesirable to restrict the 
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approach to blade models with only a few nodes at the disk-blade interface, the 
receptance method needs to be modified in order to make it even more efficient. 
This is done using the following simplifying assumption. 

It is assumed that the disk-blade interfaces undergo rigid body type translations and 
rotations The blade vibration is then determined as a combination of blade base motion 
and clamped-free blade modes As will be shown, this approach results in a formulation 
that has a relatively small number of degrees of freedom (six times the number of blades) 
and is consequently, computationally efficient. In addition, since the blade's response is 
represented in terms of its clamped-free modes its response can be quickly characterized 

with only a few modes 


Disk Equation 


If the disk substructure is subjected to a harmonic excitation, its response is given by 
the receptance matrix R, i e , 


f-d j 
< 


od 

R a,|j' 


f d l 



l>d 

_ K p,a 

pd 

R P.P. 


f d 

l f pj 


where u d , f d , and R d are the displacement, the external force, and the receptance 
associated with the disk. The subscript a denotes the group of nodes that do not interact 
with other substructures and the subscript (5 denotes the group of nodes that reside at the 
disk-blade interfaces. Since the only external forces on the disk are the interactive forces 
at the disk-blade interfaces 


Equations (1) and (2) imply that 


To formulate the disk equation in a more reduced order form, i e , a formulation with six 
degrees of freedom per interface, the foil owing two relations are introduced 
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where is a vector whose components are the six equivalent rigid body-type motions of 

'j d 

the disk at the interfaces and f£ is the resultant forces on the disk at the interfaces Q p>0 
is the geometric relation between up and uj and lj,p is the geometric relation between 
?o and fp The inverse relations of equations (4), and (5) are 


-d n d + -d 
u o=Q(3,o u p 


( 6 ) 


fd _ yd + Jd 


(7) 


d + <j + 

where Q p 0 and p are 
Their expressions are 


the generalized inverses of the matrices Q p o and I o p 


+ ^ -,d T ^d 


Qp.o - Q“m Qp.o Qp.o 


-1 


( 8 ) 


yd + _yd T f yd yd 1 

^o,p “ ^o.p I ^O.p ^o.p 


-1 


(9) 


In essence, equations (6) and (7) state that uj is the least squares fit of Up and fp is the 
non-self-equilibrated forces estimated from fjj . There are many possible choices of the 
distributed force f p The reason that only the non-self-equilibrated forces were used is 

that, because of Saint Venant's Principle, the self-equilibrated forces die off quickly away 
from the interfaces and their global effect on the system should be relatively small 
Substituting equations (6) and (7) into equation (3), results in the reduced order 
receptance formulation for the disk 


-d _ D d jd 
u o _ *^o,o * o 


( 10 ) 


where 



(11) 


d + D d 


u v u 

R o,o =Qp,o R P.P L o,P 


By applying standard modal analysis, the disk receptance R p p can be written as 

id id T 

^i.P ^.i-P (12) 


R p.p-X 


7 LJ -n 2 + 2ifto)J^ 


where «>] is the j* disk mode, and mf, , and to? are the modal mass, modal damping 
,, and natural frequency of the i* disk mode Combining equations (11) and (12), the 


ratio 


reduced order disk receptance can be written as 


id 


— h T 


R 0.0 = X 


' 1.0 ^J ,0 


d rd 


mj UJ +2iQ “j^ 


(13) 


where <j> J 0 , the j 01 reduced order disk modes, is written as 

♦io-Qp.o* 


(14) 


Blade Equation 


Consider a single blade whose base is constrained to move in a plane undergoing small, 
harmonic translations and rotations Let be the six-degree-of-freedom displacement 
vector of the blade base and, again, denote the group of nodes that do not interact with the 
disk by the subscript a and the group of nodes that reside at the interface by P The 

receptance equation for the blade can be written as 


- b r» b n b 
u a Qa,o u o 



+ Q 2 M 


b 

a t a 




(15) 


-b n b 
u p =Q 


-b 


P.o 


(16) 


where O is the excitation frequency. u b , f b , and R b are the displacement, the external 
force, and the receptance associated with the blade M b „ is the blade mass matnx 
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associated with the group-ct nodes <£.„ and Qp , 0 are geometric functions such that 
O b u b and Q» 0? describe the motion of the entire blade provided that the blade 

follows the motion of its base and does not deform The term Cl M a,ct Q a ,o u o in the 
right hand side of equation (15) describes the inertial force introduced by the blade base 
motion Rearranging equation (1 5), the motion of the group-ct nodes u„ can be expressed 
in terms of the external excitation force f £ and the blade base motion u 0 , i e. 


_b _ 0 b 
u a ” K a,a 


? b + 

1 a 


r>b 4. J 

K a,a ]Vl a,a tl 


\~b -b 
]Qa,o u o 


07 ) 


From modal analysis, the receptance R a a can be written as 


-b T 


R a.a = X 


$j.a ^j.a 


m j | (0 


b r b 


'J 




(18) 


where $ ■ is the j lh mode for a blade clamped at its base, and m • , C j . and “j are the 
modal mass, modal damping ratio, and natural frequency of the j* blade mode. 



The resultant force on the blade base, f J, is needed in order to assemble the entire 
system. To determine f J consider the blade as an object which can move and deform in 
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space (Figure 3). The forces acting on the blade are the excitation force f£ and the 
interactive force at the interface ? b . Since linear and angular momentum have to be 
conserved, one can deduce that 

^(-n 2 M b u b ) = ^f b =So.afu + fo (19) 

where X b is a geometric function which calculates the resultant force and moment about 
the blade base O. Substituting equations (16) and (17) into equation (19), f„ can be 
expressed in terms of the blade base motion u b and the external excitation f„, re., 


f b _ 7 b -b , H b f b 

f Q _ ^o,o u o +tl o,a I a 


( 20 ) 


where 


z b 0 = -£i 2 Iq M b Qo " q4 £ o,a M a,a R a,a M a,a Qa,o 


( 21 ) 


H b 0 a = - ^ So, a M l a a a , a - - o,a 


R 


-I 


( 22 ) 


Equation (20) describes the relation between the interactive force f 0 , the blade base 
motion u b , and the external excitation f b for a single blade. Clearly, similar equations 
apply for all blades In fact, equations (20), (21) and (22) describe the muttons of all of 
the blades provided the vectors represent the collections of the correspond, ng vectors of 
individual blades and the matrices become block diagonal matrices wtth the correspond, ng 
matrices of individual blades at the diagonals. 


Assembling the Substructures 


Equations (10) and (20) describe the relations between the interactive forces and the 
motion of the interfaces for the disk and the blades respectively. S.nce there are four 
unknown vectors, u d 0 , f b , u b , and f 0 b to be determined, there needs to be two add.fonal 

constraints. The first is Newton’s law of interactive forces, i e, 

fi fb (23) 

T o _ o 
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The second one is Hooke’s law which provides a constitutive relation at the interfaces, 
i.e. 



(24) 


where K 00 is the stiffness matrix associated with the attachment at the disk-blade 
interfaces 1 By solving equations (10), (20), (23), and (24) simultaneously, the four 
unknowns can be expressed in terms of the external excitation force f a , i e.. 


Uo = -(<„ + C)[l + Z* „(<„ + C)]"‘ f b 


f b = 
1 0 


+ Zo,o( R o.<, + K o 1 o)]' lH o.afa 


(25) 


(26) 


= -R 


0,0 


1 + 


7 b /»d . k -1 'll H b 

^o.o^o.o + N o,o jj n o,a 


f b 
1 a 


(27) 


f d = - 
1 o 


1 + 


Zo.o^o.o+K-'o )]" 1 


¥i b r ^ 

n o,a 1 a 


(28) 


The vibration of the blades can then be easily derived by substituting equation (25) into 
equation (17) 


Remarks on the Mathematical Approach 

There are several advantages of using the reduced order formulation (equations (17) 
and (25)) to solve the blade mistuning problem: 

1. The coefficient matrices Rq.o’ ^o,o’ anc * **o,a are calculated from the 
substructures’ modes. Since the frequency range is usually limited to a narrow 
range near a particular engine order crossing, only a few dominant modes need to 
be included in the calculation and, consequently, the computational cost is low. 


1 If the attachment is infinitely stiff, then (24) implies that the displacements in the blade and disk are 
equal, and, consequently, (24) becomes a continuity requirement. 
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2. For blade mistiming problems, the disk is treated as a cyclic symmetric structure. 
As a result, its modes can be calculated efficiently by applying cyclic symmetric 
boundary conditions on a disk segment corresponding to a single blade. 

3. The most computationally intensive part of this approach lies in computing the 
inverse of the matrix [i + Z^R^ + lC)] in equation (25). This is a square 
matrix with dimension six times the number of the blades. In many problems the 
rigidity of the disk is such that only three degrees of freedom are needed to 
describe the blade’s motion at its base, two rotations and translation normal to the 
disk. In either case the number of degrees of freedom are orders of magnitude 
smaller than for a finite element model for an entire mistuned bladed disk 

4. The reduced order method is especially computationally efficient when running 
Monte Carlo simulations of blade mistuning since the disk’s modes and the 
nominal blade’s modes need only be calculated once The coefficient matrices 
associated with the blades (z b Q , 0 and < a ) are calculated by changing the 
cantilever blade frequencies in the receptance calculation for each blade On the 
other hand, a direct simulation using finite element models would require 
completely new calculations for each mistuned stage. 


COMPARISON OF RESULTS 


A computer code named LMCC (Linear Mistuning Computer Code) has been 
developed based on the reduced order approach. A test problem was developed, Figure 4, 
and analyzed using: 1) modal summation based on a direct finite element analysis, 2) 
LMCC, and 3) the mass-spring model of Griffin and Hoosac (1984) as depicted in 
Figure 5. The response of the spring-mass model is determined using the computer code 
BLDVIB. The vibratory response of more than one hundred different resonant crossings 
have been examined of which the examples presented in this section are representative 


The geometry of the test problem was chosen so that the following requirements were 
satisfied: 

1 The disk and the nominal blade have natural frequencies that are representative of 
a turbo-pump Both the blades and disk are assumed to have material properties 

associated with a super nickel alloy. 
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2. The test problem represents a 
realistic three dimensional 
structure with low aspect ratio, 
plate-like blades. A test problem 
consisting of beam-type blades 
would not be appropriate since a 
goal is to check the validity of the 
simplifying constraint that the 
bases of the blades move as rigid 
bodies. 

3. The entire system does not have 
too many degrees of freedom As 
a result the entire bladed disk can 
be directly analyzed by finite 
elements without expending too 
much computational time or 
having numerical stability 
problems. 

There are still several parameters of 

the reduced order model that need to be 

determined before the algorithm can be 

implemented The attachment stiffness 

matrix K 00 is set to infinity to Figure 5: Spring-mass model used 

u,u 3 by Hoosac and Griffin 

represent the simplified geometry of the 

test problem. The other parameters that need to be determined are the damping ratios of 
the substructures modes. It is assumed that both the blades and the disk have the same 
damping ratio. Its value is chosen so that the resonant response predicted by LMCC for 
the tuned system is the same as that predicted by the finite element method. 

Tuned Response 

First a comparison will be made between LMCC and the finite element method based 
on the natural frequencies of a tuned bladed disk. LMCC does not readily compute 



A-11 







natural frequencies. However, if the system is lightly damped they can be inferred by the 
system’s resonant frequencies since they can be calculated using LMCC s forced 
response algorithm. Figure 6 shows a comparison of natural frequencies. The integer j 
refers to the engine order of the excitation which excites the j“> nodal diameter mode^ 
Overall, the two methods agree reasonably well However, the reduced order approach 
tends to predict slightly higher frequencies probably because of the assumption o ng. 
blade base motion which tends to over constrain the system. Note that the reduced or er 
model represents the structure’s response over a wide range of frequencies and no a 
diameters. Results from the spring-mass model are not shown since the natural 
frequencies of the tuned system are part of the input parameters to the method and not 

calculated results 



The nominal properties of the spring-mass model were calculated using the approach 
given by Griffin (1988). Initially, the spring-mass model predicted a spurious second 
mode which had to be eliminated by varying some of the parameters. Figure 1 shows 
representative response curves of the tuned system as predicted by the three different 
methods For this particular case, the bladed disk is subjected to a fourth engine or er 
excitation and the blades vibrate predominantly in their first bending mode It appears 
that the three curves agree well with each other though the response predicted by the 
reduced order model is at a slightly higher frequency. 
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Amplitudes of Blade Tips Amplitudes of Blade Tips Amplitudes of Blade Tips 



Figure 7: Tuned response curves 


Figure 8: Mistuned response curves for blades in first 

bending 
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- Mistuned Response ■ 

The finite element model is mistuned by varying the elastic modulus of the material in 
each blade. This same effect is simulated in LMCC by changing the natural frequencies of 
each blade when calculating its receptance. In the case of BLDVIB the stiffness, kjj, of 
each blade was altered to produce the same frequency distribution. Figure 8 shows 
representative response curves of the blades as predicted by the three methods In this 
case the bladed disk is subjected to a fourth engine order excitation and the blades 
predominantly vibrate in their first bending mode Clearly, the agreement of all three 
methods is quite good In general, it was observed that the reduced order model and the 
spring-mass model give comparable results when the blades respond predominately in 
their first bending modes 

Figure 9 depicts representative response curves when the blades are excited in their 
first torsional mode It is clear from the plots that the reduced order model works much 
better than the spring-mass model in representing the response of the system This case 
differs from the bending mode case in that the disk is effectively much stiffer. This 
stiffness effect is discernible from the fact that the percentage change in the frequencies of 
the four and five diameter modes of the tuned system (Figure 6) is an order of magnitude 
smaller for the torsional modes than the first bending modes. For a stiff disk, the spring- 
mass model predicts that the blades essentially respond as isolated mistuned blades on a 
rigid foundation However, both the finite element and LMCC results indicate that there 
is still significant interaction between the blades that results in complex system dynamics 
that is not captured by the simple spring-mass models 2 


Limitations of the Reduced Order Model 


Additional simulations were made with both thicker and thinner disks. In general, 
finite element and LMCC results agreed quite well except for one case. Figure 10 shows 
the response curves predicted by the three different models for a bladed disk when the 
thickness of the disk is halved In the case in question, the tuned system has a response 
which is similar to the frequency veering that is observed for the second and third modes 
when j equals seven in Figure 6. As a result, the bladed disk has two families of modes 
that are close together and that respond to the excitation. Because of the closeness of the 


2 f ' m J*! al * ons w ere also made of the response of the mistuned system in first bending with a stiffer disk. 
In this case all three methods predicted that the blades would respond as isolated blades Again, the 
spnng-mass model worked quite well for first bending. 
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frequencies, the resulting modes contain significant amounts of both disk and torsional 
blade response. It appears that neither the reduced order model nor the spring-mass 
model successfully predict the correct results. One problem with the LMCC model is 
that it does not predict the separation in the frequencies of the two families of modes 
very accurately. This problem was addressed by introducing a finite value of the blade 
attachment stiffness matrix, K^, so that the frequency difference for the tuned system’s 

response was the same as that of the finite element model. Figure 10(d) shows the 
response curves predicted by the reduced order model when the attachment stiffness has 
been adjusted. The results from the reduced order model agree somewhat better with the 
finite element results, but not as well as in the other cases examined 

— — — Execution Time Comparison 

There are two principal advantages of the reduced order model approach which lead to 
increased computational efficiency. The first is that the number of degrees of freedom 
that must be determined (three or six per blade) in the forced response analysis is 
independent of the refinement of the finite element model of the blade and disk models 
that are used to calculate its input Thus, LMCC is especially efficient when compared 
to a direct finite element solution of a bladed disk that has a large number of elements 
Secondly, the most computationally expensive part of the calculation with LMCC is in 
calculating the modes of the clamped blade and tuned disk using finite elements Since 
this needs to be done only once, LMCC becomes progressively more efficient when more 
mistuned bladed disks are simulated in a specific Monte Carlo analysis. An estimate of 
the relative run times of LMCC and a finite element analysis using modal superposition 
for a representative case indicates that LMCC would be two to three orders of magnitude 
faster for a one hundred bladed disk simulation 


Conclusions 


A new reduced order approach has been developed for analyzing the forced response 
of mistuned bladed disks. The reduced order model is based on the idea of decomposing 
the bladed disk into substructures and representing the response in terms of the degrees of 
freedom associated with their interfaces. It further assumes that the blade bases undergo 
rigid body-type motions. The blade vibration may then be expressed as a combination of 
blade base motion and cantilever blade modes. This approach results in only six 
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equivalent degrees of freedom, three rotations and three translations, at each disk -blade 
interface Because of the reduction in the number of unknowns, the reduced order model 
provides a formulation which can be solved efficiently by computers. An advantage of 
this formulation is that the reduced order model is directly calculated from finite element 

analyses of a single clamped blade and a tuned disk without any additional ad hoc 
assumptions. 


Comparisons of results from direct finite element simulations and the reduced order 
model indicate that, in general, the reduced order model agrees quite well with the finite 
element baseline Aji exception to this agreement occurred when modes from two 
different families (one predominately blade modes and one predominately disk modes) 
were close in frequency and excited simultaneously. Comparisons were also made with 
results from spnng-mass models The spring-mass models predicted good results when 
the blades were responding predominately in first bending, but did not predict the 

response when the disk was stiff and the blades were responding predominately in first 
torsion. 
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Appendix B: User Manual For LMCC 


INTRODUCTION 


LMCC (Linear Mistiming Computer Code) is a FORTRAN program designed to 
solve for the steady-state vibratory response of mistuned biaded disk assemblies It is 
the implementation of the Reduced Order Algorithm developed by Yang and Griffin 
This algorithm directly takes the results from a finite element analysis of nominal 
substructures (i.e., disk and blade) to construct a model for a mistuned biaded disk. With 
blade bases being assumed moving rigidly, it constructs a reduced order system of 
equations which can be computationally efficiently solved 

In the following descriptions of the procedure for running LMCC, it is assumed 
that the reader may refer to the example files provided with the LMCC code 


CONTROL FLOW OF LMCC 


LMCC is composed of several subprograms. Figure 1 shows the control flow of 
LMCC. The names of the subprograms under the UNIX system are shown in the 
brackets Block 1 is the interface between EAL and LMCC which rearranges the output 
data from EAL for the use of LMCC. Block 2 calculates the frequency independent 
properties of the disk and the nominal blade Block 3 helps users to specify the stiffness 
of the disk-blade attachment and generates the mistuned natural frequencies of the blades 
Block 4 utilizes the outputs from the previous blocks to calculate the forced response of 
the mistuned biaded disk and generate preliminary results from which the vibration of the 
disk and the blade can be calculated Block 5 extracts the information from preliminary 
results to calculate and print out the nodal amplitudes of the blades. 

An important rule to run LMCC is that the subprograms in the same block are 
mutually independent but are dependent of the results from the previous blocks: 
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Example 1. If a user wants to change the attachment stiffness of the disk -blade 
interface but keep analyzing the same set of blades after a certain run, he/she can 
just go back to Block 3 to run the subprogram [interface] and run [assembly] in 
Block 4 and [amplitude] in Block 5, consequently; he/she does not have to run 
[ranbld] in Block 3 or any subprogram in Blocks 1 and 2. 

Example 2 After a certain run, if a user decides that he/she wants to see the 
amplitudes of certain blade nodes which were not printed out previously, he/she 
only needs to rerun the subprogram [amplitude]. 


Example 3. To run a Monte Carlo simulation for the mistuned analysis of a particular 
nominal design, one can run LMCC according to the following steps 


Step 1 

Run EAL for the disk and the nominal blade 

Step 2. 

Run [eal_disk] and [eal blade] 

Step 3 

Run [disk] and [blade] 

Step 4. 

Run [interface]. 

Step 5. 

Run [ranbld] to generate mistuned blade frequencies Giving different 
seeds to the random number generator will produce different sets of 
blades 

Step 6 

Run [assembly]. 

Step 7. 

Run [amplitude] 

Step 8. 

If necessary, rename the output file of [amplitude] to save the result of 
a particular stage The default output file name is 
Imcc. out. amplitude ” If this file is not renamed, the next run of 
[amplitude] will overwrite the file. 

Step 9. 

Change the seed of [ranbld] and go back to Step 5. 
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Read in disk data, 
[eal disk] 


Read in blade data, 
[eal blade] 


Calculate frequency 
independent terms for 
the disk. 

[disk] 


Calculate frequency 
independent terms for 
the nominal blade, 
[blade] 


Specify attachment 
stiffness between the 
disk and the blade, 
[interface] 


Generate mistuned 
blade frequencies, 
[raribld] 


Calculate forced response of the mistuned 
bladed disk over a specified frequency range 
and generate preliminary result. 

[assembly] 


Print out nodal arrplitudes of the blades, 
[amplitude] 


Figure 1 Control Flow of LMCC 
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INSTALLATION OF LMCC 


To install LMCC under UNIX systems, one can do the following: 


1. Store all of the LMCC source codes in a directory and cd to this directory. 

2. Key in make to create the executables automatically. You can also build LMCC 
manually by keying in the following commands: 

f 77 -o eal_disk eal_disk.f 
f 77 -o eal_blade eal_blade.f 

f 77 -o disk disk_0 . f disk_l.f disk_2.f dmatrix.f zmatrix.f 

f 77 -o blade blade_0.f blade_l.f blade_2.f dmatrix.f zmatrix.f 

f 77 -o interface interface. f 

f 77 -o ranbld ranbld_0.f ranbld_l.f 

f77 -o assembly -O assembly. f dmatrix.f zmatrix.f 

f 77 -o amplitude -0 amplitude. f 


Notice : Since LMCC is quite computationally intensive, users are encouraged 

to compile the optimized executables. However, due to problems from 
some vendors, the optimization should be used with care. For 
example, on Silicon Graphic workstations running IRIX 5.2, using 
MIPS FORTRAN 77 to compile ranbld with optimization can cause 
the executable to work improperly 


3. Correctly set the PATH environment so that the executables of LMCC can be 
found by the system. 


INPUT FILES FOR EAL 

Before running LMCC, the finite element analysis by EAL has to be done first 
The format of the EAL input files for the disk and the (nominal) blade are as follows: 

Input File for the Disk 

The input file for the disk must contain the following 
statements (in addition to the standard input data, for instance, 
geometry, element type, etc.) This version of LMCC requires that a 
modal analysis of the entire disk is made, that is, not a cyclic 
symmetric segment. 
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*XQT EIG 


$ END OF EIG 
*XQT DCU 
PRINT 1 JLOC 
PRINT 1 VI BR MODE 


$ END OF DCU 


Input File for Ih& (Nominal) Blade 

The input file for the blade must contain the following 
statements (in addition to the standard input data, for instance, 
geometry, element type, etc.) 


*XQT M 


$ 


END OF M 


*XQT EIG 
RESET M=CEM 


$ END OF EIG 
*XQT PS 
*ONLINE=l 
*FORMAT=2 
CEM 

$ END OF PS 
*XQT DCU 
PRINT I JLOC 
PRINT 1 VIBR MODE 


$ END OF DCU 
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Notice 


, The execution sequence shown above should be strictly followed. Othenv.se 
LMCC will indicate there is an error. The general rule is that 'he n», oral 
frequencies have to be printed out firs, (by EIG). then CEM matnx (by PS), then 

JLOC (by DCU), and finally VBR MODE (by DCU). 

2 LMCC needs the complete mass matrix of the nominal blade to calculate the 
forced response of the system. Since EAL can no. print out the moments of 
inertia of the nodes while using lumped mass matrix (DEM), LMCC uses 
consistent mass matrix (CEM). Therefore, when runnmg EAL on the blade, one 

must use CEM for the analysis. 

3 The join, locations (JLOC) and the natural modes (VIBR MODE) have to be pu, 
in Cartesian coordinates LMCC will not work properly ,f they are specified in 

polar coordinates 



To run the subprograms of LMCC, the associated input files have to be ready 
beforehand. There are eight input files corresponding to each of the 
indicated in Figure 1 The input files are named as [ subprogram namej.m ^ 

the input file name for subprogram eat disk is ealjhsk.m e inpu 
information which are essential to the forced response analysts and wrll be read by 
subprograms The formats of the input files are as follows: 


Example File 1: Input File for Subprogram eal.disk - eal dWc.fn 

The subprogram eatjisk functions as a interface between EAL and LMCC w 
reads the EAL ouput for the disk and reformat it for the input of LMCC 


1.0 

0.0 

0.0 

0.866 

0.5 

0.0 


0.0 

0.0 

1.0 

0.0 

0.0 

1.0 

-0.5 

0.0 

0.866 

0.0 

0.0 

1.0 


•.record 1 
: record 2 
: record 3 

: record 4 
: record 5 
: record 6 
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NDI: record 7 

3 



: record 8 

481 



: record 9 

482 

483 



: record 10 (7+NDI) 

0.0 

5.0 

0.0 

: record 11 (7+NDI+l) 

10 



: record 12 (7+NDI+2) 




•.record 13 (7+NDI+3) 

45 



:record 14 (7+NDI+4) 

12 

3 

4 

5 

:record 15 (7+NDI+5) 


F.xplanatip n of Records 

(1 to 3) . r , 

Matrix of direction cosines which transfers the disk coordinates to the first 

blade coordinates 


(4 to 6) ,. . 

Matrix of direction cosines which transfers the ith blade coordinates to e 
(i+l)th blade coordinates. 

( 7 ) . ^ 

NDI: Number of disk nodes per disk-blade interface. 

(8 to 7+NDI) . 

Nodal numbers of disk nodes residing at the first disk-blade interface 

Position of the reference point O (in disk coordinates) assoc, ated wtth the first 
blade 


(7+NDI+2) ... 

Nodal number difference between corresponding disk nodes residing a 

neighboring interfaces 

(7+NDI+3) 

Number of disk modes to be stored in the library (for future use). 

(7+ND+4) 

Number of blades. 

(7+ND+5) 

Degrees of freedom of VIBR MODE printed out by DCU 


File eal_blade.in - Input File for Subprogram eaU>I«fe 


10 


: record 1 
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2 NBI : record 2 

1 : record 3 

~ : record 4 

2 : record 5 (2+NBI) 

0.0 5.0 0.0 : record 6 (2+NBI+l) 

245 : record 7 (2+NBI+2) 

Explanation of Records 
(1) 

Number of blade modes to be stored in the library (for future use) 

( 2 ) 

NBI: Number of blade nodes per disk-blade interface 
(3 to 2+NBI) 

Nodal numbers of blade nodes residing at disk -blade interface. 
(2+NBI+l) 

Position of the reference point O (in blade coordinates). 
(2+NBI+2) 

Degrees of freedom of VIBR MODE printed out by DCU. 


File disk An — Input File for Subprogram disk 


10 

2 

3 

8 

9 

13 

14 
21 
22 
23 
25 


NDM: record 1 
: record 2 
: record 3 
: record 4 
: record 5 
: record 6 
: record 7 
: record 8 
: record 9 
•.record 10 
: record 11 (1+NDM) 


Explanatio n of Records 

0 ) 

NDM: Number of disk modes to be included in forced response calculation 
(2 to 1+NDM) 

Mode numbers of disk modes to be included in calculation 
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Pile blade.in — Input File for Subprogram blade 


3 

1 

3 

4 


NBM: record 1 
: record 2 
: record 3 
: record 4 (1+NBM) 


2 

8 3 (1.0, 0.0) 

9 5 (1.0, 0.0) 


NF: record 5 (1+NBM+l) 

: record 6 (l+NBM+2) 

: record 7 (1+NBM+l +NF) 


Explanation of Records 

( 1 ) 

NBM: Number of blade modes to be included in forced response calculation 
(2 to 1+NBM) 

Mode numbers of blade modes to be included in calculation 
(1+NBM+l) 

NF: Number of non-zero components of excitation forces to be specified 
(l+NBM+2 to 1+NBM+l+NF) 

Non-zero components (in blade coordinates) of excitation forces In this 
example, record 6 specifies that force component 3 on blade node 8 is (1 .0,0.0) 
and record 7 specifies that force component 5 on blade node 9 is (1 .0,0.0). 


Notice that force components have to be specified by complex numbers to 
take into account the possible phase difference between force components 


File interjace.in — Input File for Subprogram interface 


6 

1 

1 

(1.00E55, 

0.0) 

NK: record 1 
: record 2 

2 

2 

(1.00E55, 

0.0) 

: record 3 

3 

3 

(1.00E08, 

0.0) 

: record 4 

4 

4 

(8.70E05, 

0.0) 

: record 5 

5 

5 

(7.20E08, 

0.0) 

: record 6 

6 

6 

(1.00E55, 

0.0) 

: record 7 (1+NK) 


Explanation of Records 

( 1 ) 

NK: Number of non-zero components (in blade coordinates) to be specified 
for the disk-blade attachment stiffness matrix KoO- (Koo is a 6x6 complex 
matrix.) 
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(2 to 1+NK) 

Non-zero components of Koo In this example, 
record 2: stiffness in x direction is infinite, 
record 3. stiffness in y direction is infinite, 
record 4. stiffness in z direction is (1.00E08,0.0). 
record 5. torsional stiffness in x direction is (8.70E05,0.0) 
record 6: torsional stiffness in y direction is (7.20E08,0.0). 
record 7: torsional stiffness in z direction is infinite 


Nolkfi that stiffness components have to be specified by complex numbers to 
take into account the possible damping mechanism. 


File ranbld.in — Input File for Subprogram ranbld 


3 

1 

3 

50.0 

874377 


: record 1 
: record 2 
: record 3 
: record 4 
: record 5 


Explanation of Records 

( 1 ) 

Indicator for blade natural frequency disturbance. 

1 = identical blades. 

2 = sinusoidal distribution. 

3 = gaussian distribution 

( 2 ) 

Blade mode to be disturbed. 

( 3 ) 

Indicator for the frequency disturbance on other modes. 

1 = no disturbance 

2 = same magnitude 

3 = same percentage 


Amplitude (sinusoidal distribution) or standard deviation (gaussian 
distribution) of the disturbance in Hz 
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Number of nodal diameters of sinusoidal distribution or seed for the gaussian 
random number generator 


Mle assembly. in - Input File for Subprogram assembly 


1 


: record 1 


4 

4450.0 

4950.0 100 

: record 2 
: record 3 


0.0025 

2 

13 

14 

0.0020 

0.0015 

•.record 4 
ND: record 
: record 6 
: record 7 

5 

(5+ND) 

0.0027 

2 

3 

4 

0.0022 

0.0024 

•.record 8 
NB: record 
•.record 10 
: record 11 

(5+ND+l) 

9 (5+ND+2) 
(5+ND+3) 
(5+ND+2+NB) 


F.xplanatio n of Records 

0 ) 

Flag for error estimate on LU algorithm. 

1 = perform 

0 = do not perform 

( 2 ) 

Excitation engine order 

Excitation frequency bounds (in Hz) and number of the intervals of the 
calculation In this example, the lower bound is 4450 Hz, the upper bound is 
4950 Hz, and the number of intervals is 100. The forced response will be 
calculated for every 4 Hz from 4450 Hz to 4950 Hz 

(4) 

Damping ratio for disk modes in general 

( } ND: Number of disk modes with damping ratios to be particularly specified. 
(6 to 5+ND) 

Mode numbers and damping ratios of the particular disk modes. In this 
example, record 6 specifies that the damping ratio for disk mode 13 is 0.0020 
and record 7 specifies that the damping ratio for disk mode 14 is 0.0015. 

(5+ND+l) 
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Damping ratio for blade modes in general. 
(5+ND+2) 

NB: Number of blade modes with damping 


ratios to be particularly specified 


(5+ND+3 to 5+ND+2+NB) 

Mode numbers and damping ratios of the particular blade modes. In this 
example, record 10 specifies that the damping ratio for blade mode 3 is 0 0022 
and record 1 1 specifies that the damping ratio for blade mode 4 is 0 0024. 


File amplitude.in - Input File for Subprogram amplitude 

NN: record 1 

\ : record 2 

^ : record 3 (1+NN) 

8 

NC: record 4 (1+NN+l) 

^ : record 5 (l+NN+2) 

3 -.record 6 (1+NN+l+NC) 


Explanatio n of Records 

( 1} 

NN: Number of blade nodes to be printed 
(2 to 1+NN) 

Blade nodes to be primed In this example, amplitudes of nodes 3 and 8 wtll 

be printed 
(1+NN+l) 

NC: Number of displacement components to be printed for each no e 

specified above. 

(l+NN+2 to 1 +NN+l+NC) 

Displacement components (in blade coordinates) to be printed. In this 
example, displacement components 3 and 5 will be printed. 



After successfully running subprogram amplitude , the forced response of the 
bladed disk is stored in file Imcc, out, amplitude This file has documents in the beginning 
which explain the format of the data printed. The format is designed to be read by general 
plotting routines. An example of the file is as follows: 
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# Data printed are: 


# (column 1) Excitation frequency [Hz] 


# 

(column 

2) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

1 

# 

(column 

3) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

2 

# 

(column 

4) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

3 

# 

(column 

5) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

4 

# 

(column 

6) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

5 

# 

(column 

7) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

6 

# 

(column 

8) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

7 

# 

(column 

9) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

8 

# 

(column 

10) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

9 

# 

(column 

11) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

10 

# 

(column 

12) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

11 

# 

(column 

13) 

Amplitude 

of 

component 

3 

of 

node 

4 

of 

blade 

12 

# 

(column 

14) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

1 

# 

(column 

15) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

2 

# 

(column 

16) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

3 

# 

(column 

17) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

4 

# 

(column 

18) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

5 

# 

(column 

19) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

6 

# 

(column 

20) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

7 

# 

(column 

21) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

8 

# 

(column 

22) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

9 

# 

(column 

23) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

10 

# 

(column 

24) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

11 

# 

(column 

25) 

Amplitude 

of 

component 

3 

of 

node 

8 

of 

blade 

12 

# 

(column 

26) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

1 

# 

(column 

27) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

2 

# 

(column 

28) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

3 

# 

(column 

29) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

4 

# 

(column 

30) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

5 

# 

(column 

31) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

6 

# 

(column 

32) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

7 

# 

(column 

33) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

8 

# 

(column 

34) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

9 

# 

(column 

35) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

10 

# 

(column 

36) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

11 

# 

(column 

37) 

Amplitude 

of 

component 

5 

of 

node 

4 

of 

blade 

12 

# 

(column 

38) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

1 

# 

(column 

39) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

2 

# 

(column 

40) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

3 

# 

(column 

41) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

4 

# 

(column 

42) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

5 

# 

(column 

43) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

6 

# 

(column 

44) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

7 

# 

(column 

45) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

8 

# 

(column 

46) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

9 

# 

(column 

47) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

10 

# 

(column 

48) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 

blade 

11 

# 

(column 

49) 

Amplitude 

of 

component 

5 

of 

node 

8 

of 
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additional consideration and potential problems 


While using EAL as the preprocessor of LMCC, we found that EAL occasionally 
generates disk modes which are not mutually orthogonal. The problem was that modes of 
the same number of nodal circles and same number of nodal lines should have the same 
mode shape except that they are out of phase. EAL sometimes generates this kind of 
modes with correct mode shapes which are not completely out of phase This kind of 
error causes the forced response calculation so that LMCC will give meaningless results 
Typically when this error occurs, it is indicated by the fact that LMCC predicts 
unusually high forced response over part of the frequency range. 
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Appendix C: Theoretical Manual for LMCC 



This appendix contains the theory of the reduced order algorithm upon which the Linear 
Mistuning Computer Code (LMCC) is based Chapter 1 is the introduction and motivation of 
this research project. Chapter 2 and Chapter 3 state the theory of the reduced order 
algorithm. Chapter 4 gives some benchmark test cases and the comparison between t 
reduced order model and the traditional spring-mass model. 


AND MOTIVATION 


In the use of turbomachineiy in modem aircraft, blades have failed from high cycle fatigue 
because of excessive vibration Based on experiments and earlier analytical investigations 
[1-5] it is realized that, because of the rotational periodicity of its geometry, a bladed disk 
usually has natural frequencies that are clustered in narrow ranges. When the natura 
frequencies of a system are close together, slight variations in the system’s structural 
properties can cause large changes in its modes, and, consequently, its dynamic response. The 
sensitivity of a bladed disk’s dynamic response to small variations in the frequencies o e 
blades is referred to in the literature as the “blade mistuning problem” and has been studied 
extensively It is important to understand mistuning since it can result in large blade to blade 
variations and the high response blades can fail from high cycle fatigue. For modem aircraft, 
such as a space shuttle, the high vibratory response of bladed disks is of special concern 
because of the severe operating conditions to which they are subjected 


Every year, a great number of experiments are conducted by industry to ensure that the 
ki,u». are safely designed These experiments are generally expensive and time-consuming In 
order to reduce the time and cost associated with these experiments, turbine manufacturers 
replace part of the experiments by numerical simulations. From a modeling point of view, 
traditional simulations can be classified into several categories One of them is to directly use 
finite element methods The direct finite element approach is straight forward and can 
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accurately model the characteristics of a .^lemenr mriel for a 

with enough elements. However, since a suite ^ s rf freed om, the associated 

typical bladed disk can have more tha , , ^ Qf , si „ gle suge can be as 

computational time for denvtng “ m istuning ana i ysis usi „g Monte Carlo 

high as 1,000 CPU minutes on CRAY ^ ^ ^ be prohlb i,ively high. (For 

Simulation is desired, the simulation.) In addition, the extremely large 

r r::;= r — - — 

instabilities make the direct finite element approach unfeas. e. 

There have been attempts to develop i^'ch^e W^s in modeling 

Reference may be made to a feature y rf , he ^mpu to develop structurally 

blades for vibration analyses represent the disk and beam elements 

accurate models ™d disks use ^ conflgur a,ions for which the 

to represent the blades [3,1 1-13,2 1 ^ eve „ solid elements are often 

beam representation may be adequate, p a , wuh recent improvements in finite 

needed to represent modem, low aspec ra ■ ^ fee mode led accurately but, 

element methods and computer spee , a w ^ ^ of ^mg th ese programs would be 
as previously mentioned, tt is recogmz ^ ^ , 0 , he attem p,s by researchers 

prohibitively high for use dunng the desig p m fini ,e element analyses of individual 

to use substructuring techniques in which results ^ ^ ^ ^ bladed disk [14 , 151 

substructures are used to obtain the ynamic^ symmetric boun dary conditions into 

There have also been attempts to intro u ^ m modes ^ be effic iently solved 

"rr — Ir^uencies being sepamted, the numerical stabili^ is 
improved [7,16,17] 

. u with cvclic symmetric boundary conditions can be easily 
While finite element models y been done in mistuning utilizes 

solved for tuned bladed disks, much o t e wor reduce the number 

spring-mass models to represent the bladed disk IW whi , e they 

systems. Little work has been published which verifies their a«u 
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The goal of this research is to develop a methodology that will directly take the results 
from a finite element analysis of a tuned system and construct a computationally efficient 
mistuned model with a reduced number of degrees of freedom. The intent is that the approach 
will display the structural accuracy of the finite element model and computational efficiency 
comparable to a spring-mass model 


Chapter 2: Receptance Method And Reduced Order 
Modeling 


— — 2.1 Receptance Method 

In structural dynamics, one widely used analytical method is to subdivide a complex 
structure into several substructures, analyze the behavior of the individual substructures, and 
then use the receptance method [6] to reassemble the structure. Receptance method is based 
on the observation that, for every substructure, the dynamic response is determined by how 
the substructure interacts with its environment. If the substructure interacts with its 
environment only at a limited number of interfaces, it is possible to express the degrees of 
freedom of the entire substructure in terms of the degrees of freedom of these interfaces. The 
benefit of this method is that when several substructures interact with each other, we only 
need to solve for the degrees of freedom of the interfaces. Once the degrees of freedom of 
these interfaces are determined the response of all substructures and, therefore, the whole 
structure is known 


To apply the receptance method to a bladed disk, let us consider the disk and all of the 
blades as two separate substructures We can then write the receptance equations for the disk 
and the blades 


- d 
u 



( 2 . 1 ) 


- b 
U 



( 2 . 2 ) 


0 ° H 

where u , f , and R° are the displacement, the external force, the receptance associated with 

b ** b u 

the disk, and u , f , and R are the displacement, the external force, and the receptance 
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write equations (2 1) and (2 2) in components, 


can 




R d R° 

R a, a R a,P 

|«1 

i-d| 
\ u p 1 


pd j»d 

l K p,a K P.PJ 

Up 

- b 1 
U a 

' 

l. 

T»b |>b 

R a,a" a,P 

I - b \ 

M 

— b 

\ up 


p b D b 

[R P,a K P.Pj 

U 


(2.3) 


(24) 


Since, in general, the only external forces on 
interfaces, we assume that 


the disk are the interactive forces at the disk-blade 


fa=0 


(25) 


requiring the displacements to be continuous and the interactive to 
Interactive Force Law, 


— d ■* b 

U p = U p 


-d ■* b 
f p = -f P 


( 2 . 6 ) 

(27) 


we can solve for the forces at the interfaces using equations (2.3) and (2.4) 

f P = -(rVp + |>,J ^ 


and 


( 2 . 8 ) 
(29) 

Since the only unknowns on the right hand sides of equations (2.3) and (2.4) are f p a" d f P • 
the response of the entire bladed disk is known. 


fp =(R b p.P +R Vp) R V Ja 


■* b 
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The benefit of using equations (2 8) and (2 9) is that we only need to solve the problem in 
the subspace-p without the necessity of either calculating the inverse of the dynamic 
impedance of the whole bladed disk or finding the modes for the entire system. However, 
finding the receptances for the disk and the blades beforehand still can be computationally 
intensive. A mathematical expression for the receptances of the disk and the blades are 

R d =<K d - Q 2 M d + in C * 1 ) -1 (2 ]0) 

R b =(K b -ft 2 M b + ifiC*)"' (2I)) 

where K d , M d , and are the stiffness, mass, and damping matrices of the disk, and K b , 
M b , and are the stiffness, mass, and damping matrices of the blades £2 is the excitation 
frequency. R b is a block-diagonal matrix if there is no structural links between the blades 


An alternative way of calculating these receptances is to use the substructures’ modes 
That is, by applying standard modal analysis, the receptances can be written in the following 
form 


Rd 4> j 4>j 

j m d j (co d j - Q 2 + 2 i Q toj £ d 

-b-b T 

Rb _y 0 j <t>j 

j b | b 2 2 ~ „ b r b 

J m j (co j -Q +2i£2cOjQj 

d d d d 

where <f> j , m :, L ;, and co ; are the mode shape, modal mass, modal damping ratio, and natural 

J J J "* b b b b 

frequency of the j-th disk mode, and <|>j , m :, £j, and co j are the mode shape, modal mass, 
modal damping ratio, and natural frequency of the j-th blade mode, respectively. This method 
is in general better that equations (2.10) and (2.11) since, with the substructures’ modes 
derived in advance, using equations (2 1 2) and (2. 1 3) to calculate the receptances needs merely 
algebraic multiplications and additions, whereas equations (2.10) and (2.11) involve the 
calculation of matrix inverses for every excitation frequency. The other benefit of equations 
(2.12) and (2.13) is that not all of the substructures’ modes are needed in the calculation. 
Since the modes with natural frequencies close to the excitation frequency will dominate the 
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response, the other modes can be eliminated from the calculation without a significant loss of 
accuracy. 


However, there are still several shortcomings which prevent us from using the receptance 
formulation (the combination of equations (2.8), (2.9), (2.12), and (2.13)) to solve our 
problem: 

1. The substructures’ modes used in equations (2.12) and (2.13) have to be free at the 
disk -blade interfaces, or, in other words, they have to be admissible. To satisfy this 
requirement, it is fine for us to use the disk modes which are clamped at the inner edge 
and free at the outer edge, but it is undesirable to use the free-free blade modes. 
Because the blades generally vibrate close to the clamped-free condition, we need a 
large number of the free-free modes to achieve a good representation of blade vibration 
Another more realistic reason for not using free-free blade modes is that most of the 
analyses done on the blades in the industry, either numerically or experimentally, use 
cantilever blades There is much more information available on the clamped-free blades 
than the free-free blades 

2. In equations (2.8) and (2.9), we need to find the inverse of the matrix (r + R p>p ) 
This matrix has a dimension of the group-P nodes which is much smaller than the 
dimension of the whole bladed disk, but, still, there can be too many degrees of 
freedom to be solved for Since we do not want to restrict our analysis to use a 
particular type of blade models, for example, beam-type blade model with single-node 
contact at the disk -blade interface, we need to modify the receptance method to have 
the capability to handle multiple-point disk-blade attachments. 


— — i 2.2 Reduced Order Modeling Based On Rigid Blade Base 

Assumption — 

One of the reasons that the direct receptance method is not particularly attractive is that 
there are still too many degrees of freedom. In this section, we would like to propose a 
reduced order method which is capable of handling multiple-point disk-blade attachments and, 
yet, results in a formulation having fewer degrees of freedom. 
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For a system of substructures having multiple point contact with each other, . ,s tn 
general not possible to reduce the number of degrees of freedom ofits receptance formu, attorn 
However for a bladed disk, the defotmation of the disk-blade interfaces ,s tn general 
"™f,cam in comparison with the deformation of the blades. Therefore, except for the , cas 
when the high nodal diameter disk modes strongly participate, each o t e a e ases 
reasonably* approximated by six degrees of freedom, namely, three translations and three 
rotations- With this assumption, the number of degrees of freedom o * ' ” ° 

formulation will be six times of the number of blades whtch ,s much smaller than 

direct receptance formulation 



Figure 3 1 shows a schematic of a bladed disk assembly. By subdividing the structure 
along the dashed line as indicated in the figure, we have two types of : subt ttrucntres, mtuney^ 

the disk (Figure 3.2) and the blades (Figure 3.3). Observing Figures an ., 
matte Imion of the disk is a function of the interactive forces a, the disk-blade interfaces 
in addition, the dynamic response of the blades is also determined by these — « «« 
provided the external excitation on the airfoils are taken tnto account. As result, 
vibration of the entire structure can be expressed in terms of the interactive forces 



Figure 3-1 : A schematic graph for a bladed 
disk assembly 
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The goal of this chapter is to present the 
application of the described reduced order 
approach to bladed disk assemblies. In 
Sections 3.1 and 3 2, we will present the 
mathematical formulations which express the 
vibration of the disk and the blades in terms of 
the reduced order degrees of freedom of the 
disk-blade interfaces. In Section 3.3, we will 
derive the solutions for the whole structure by 
utilizing the results of Sections 3.1 and 3.2 



Figure 3-3 : Substructure blades 


3.1 Formulation for the Disk 


3.1.1 Disk receptance 

For the disk substructure which is clamped at the inner edge and free at the outer edge, the 
equation of motion is 


M d x d +C 1 x* +K d x d = f d e int 

where M d , C 3 , and K d are the mass, damping, and stiffness matrices respectively, x is the 

complex displacement vector whose real part is the physical displacement of the disk, A is the 

** d 

excitation frequency, and f e is the complex forcing vector whose real part is the 
physical force acting on the disk. When the bladed disk is operated at a stationery speed, the 
steady-state response of the disk can be written as 



- d . 

where u is the complex amplitude vector of the response By substituting (3.2) into (3.1) 
and eliminating time dependent terms, we obtain 


- d 
u 



(3.3) 
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where R d is the receptance matrix of the disk with the expression 

R d =(K d — a 2 M d + iOC* 1 ) (3.4) 

By applying standard modal analysis technique and approximating the structural damping by 
modal damping, the disk receptance can be further written as 


R d =X 


-* d d 

<> j <i>j 


i m d (a> d 2 - fi 2 + 2 i ft ti)j C d j) 


(3-5) 


whereof is the j-th disk mode, and m d , and m d are the modal mass, modal damping ratio, 
and natural frequency of the j-th disk mode, respecuvely 

To calculate the receptance R d , equatton (3.5) is preferable than equation (3.4) ftom the 
computational point of view since equation (3.5) inheres several mce propen, es. 

1 Equation (3.4) involves the calculation of a matrix inverse which is in general 

computationally expensive. In the case of near or a. the natural frequency of the d^ 

the use of equation (3.4) is panicularly undesirable smce die dynamic impedance 
/ ^ j .... rnntrarv. eauatlOH (3.3) 


on (3 4) is particularly undesirable since me uy— —r 

( K d _o 2 M d + in C 5 ) tends to be ill-conditioned On the contrary, equation (3.5) 
needs metely algebraic multiplications and additions I, is therefore computationally 

more efficient and reliable. 

2 For a particular frequency range of interest, only a limited number of the disk modes 
dominate the response and the other modes participate in an insignificant mam, 
eliminating the insignificant disk modes ftom the calculanon equation a t, 
remains an accurate representation of the disk reliance and yet 1 .as the "ageo 
higher computational efficiency. On the other hand, using equation (3 4) to c^culMe 
the disk receptance would require an intensive computational effort independe 
frequency range of interest, 

3 The disk receptance for the teduced degrees of freedom of the disk interfaces can be 
easily derived using equation (3.5) and tesults in a computationally efficient appro * . 

However, using equation (3.4) for the same purpose does not have such a bene , 
will see the reason for this in the next section. 
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3.1.2 Disk receptance at the disk-blade Interfaces 

As was mentioned in Chapter 2, the response of the disk is a function of the interactive 
forces that act between the disk and the blades. To derive the disk receptance associated witjj 
the disk-blade interfaces, we need to write the force vector f and the mode shape vector <f) j 
in components, 


- d 

♦ j 


- d 


4* j,a 



f 


d 



- d 
fa 


- d 
f P 


(3.6) 


(3.7) 


where the subscript “P” indicates the disk -blade interfaces and “cc indicates the other parts 
of the disk. Since the only forces acting on the disk are the interactive forces between the disk 
and the blades, we have 



(3.8) 


Combining equations (3.3), (3.5), (3 .6), (3.7), and (3.8) we obtain 


- d 
U 





Qj 4>j,p 

- Q 2 + 2 i Q 0) d j C d j) 


- d 

fp 


(3.9) 


In particular, if we are interested in the displacement of the interfaces, equation (3.9) implies 


— d d "* d 

u P = R p,p f P (3.10) 

where u p is the displacement of the disk at the interfaces. R pp is the disk receptance 
associated with the interfaces and is given by. 
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R Vi! = I 


• d -d T 
<t> j,p j.P 


> m d j (<o d j 2 -n 2 + 2iO“jC'j) 


(3 H) 


The key results of this section is the reduced disk receptance equation 

*. relation between the forces and the displacentents of the nodes assoc, a, ed w„h the d.sk 
interface degrees of freedom 

3.1.3 Deriving sir enuiuaient degress el iisalam lat sash dlsk-ft'atfe l "** rla '” 

I, was discussed in Chapter 2 that the disk receptance associated w,th the interfaces st 
has too many degrees of freedom to be optimally efficient. The goal of tins section is t 
derive a reduced ordered disk receptance which has significantly fewer degrees of freedom^ 
The basic idea is to find six equivalent degrees of freedom which will approximately repre 
"Ich disk -blade interface Assume that, fo, the j-th 

and rotation about a reference point O is approximately given by v . and . «p V 
(Here, we omit the subscript T temporarily for the simplicity of wnhng ) f we deno e ^ 
translation of the k-.h disk node of this interface as v k and the rotation of the node k, 
then, for small amplitude vibration, we can wnte 


-d - d ;* d - d 

Vk~ v o + ®o x Pk 


(3.12) 


- d - d 
0 k = 0 o 


(3.13) 


- d 


where p k is the position vector 


from point 0 to the k-th node. If we define the vectors 


- d 
U o = 


- d 
U k = 


* d 

1 e o 

- d 
V k 


(3 14) 


I 


- d 

6k 


(3.15) 


k = 1, 2, .... DN1, where DN1 is the number of disk nodes at the disk-blade interface, then 
equations (3.12) and (3.13) can be cast in matrix form 
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where 


Q d k, 0 = 


i w 


0 


k,o 


(3.17) 


W' 


k,o 


o 

— d 
“IP k 
d 


h A 

h d l 

lp klz 

“IP k ! y 

0 

fel 

-in 

0 


(3.18) 


and I is the 3 by 3 identity matrix. 


- d 


Now, since u p is the displacement vector for all of the disk nodes at the disk-blade 
interface, we have 


- d 1 
up = 


** d 1 -* d 1 ■* d 

U 1 ...U k ...U dni 


(3.19) 


Then equations (3.16) and (3.19) imply 

- d f-wd d 

U P " Qp,o u o 


(3.20) 


where 


Qp.o = 


Q l,o -Q k,o -Q dNI.0. 


(3.21) 


Without dropping the subscript “j” for the j-th disk-blade interface, equation (3.20) would be 
written as 


(3.22) 


(u p) j = (Q d p 0 )j (u o)j 

Then, for all of the interfaces, equation (3.22) implies 



(3.23) 


- d 


u p = Q 


P.o 


d 

o 


where, with the number of blades being denoted by NB, 


- d 


T 


Up 




T 

NB 



(3.24) 

(3.25) 


- d 1 
U o = 


- d 1 
u of 1 


T U 


U of k ...\U o/NB 


(3.26) 


Notice here that Q p >0 is a block-diagonal matrix with diagonal blocks (Q d p J k (k = 1 , ..., NB) 

Equation (3.23) gives an approximate relation between the reduced order displacements u 0 

d d 

and the displacements of the interfacial disk nodes up. u 0 is a vector with dimension 6«NB 

-* d (j 

and u p has a dimension of b^DNBNB Qp 0 is a rectangular matrix with dimension 
6*DNI»NB by 6»NB 


We also need to find the resultant forces and moments for the each interface. For the j-th 
disk interface, the resultant force and moment about the reference point O are 


- d 
n o 



k 


(3.27) 


- d _ - d — d — d 

'to = 2. t k + PkXTt k 

k 


(3.28) 


where n 0 is the resultant force acting on the disk at the disk interface and 1 0 is the resultant 

~ * d 

moment about the reference point 0 acting on^tlje disk at the interface. K k is the force acting 
on the k-th disk node of the j-th interface and t k is the moment acting on the k-the disk node 
of the j-th interface. Notice here that the subscript “j” is temporarily omitted for simplicity in 
writing. Equations (3.27) and (3.28) can be rewritten in matrix form by using the procedure 
presented earlier, 



0 

I 


- d 

n k 

- d 
X k 


(3.29) 
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where the skew-symmetric matrix W k 0 has been defined in equation (3.1 8). If we define the 
vector 


f d J *0 1 

!;:f 


(3.30) 


and recognize that 


-d T 
f P = 


-d T -d T 
7t 1 X 1 


T.hT _ 


T _ 


.7t k T k --Tt DNI T DNI 

then equation (3.30) can be written as 
- d d T - d 

f o= Qp , 0 fp 


(3.31) 


(3.33) 


where Qp c has been previously defined in equation (3.21). Now, using the subscript “j” 
which was dropped previously, equation (3 33) would be in fact written as 




Therefore, for all of the disk-blade interfaces, equation (3.34) implies 
*• d d T -» d 

f o = Qp , 0 f P 


where, with NB denoting the number of blades, 


-d T 

f o — 

-d T 

4 = 


(3.34) 


(3.35) 


h “1 T 

b d l T fc d l 

T 

\f of 1 — 

f oj k 4f o, 

NB. 

(; d ) T i 

M T /; d l 

i T 

]f p] i ...1 

i f p)k -if 

'nb. 


(3.36) 
(3 37) 


and Q p c was defined in equation (3.25) Equation (3.35) governs the relation between the 

resultant interfacial forces, f 0 , and the individual nodal forces, f of the interfacial nodes for 

■* d ■* d K 

every disk interface, f 0 and f p are vectors with dimensions 6*NB and 6*DNI*NB 
respectively. 
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3.1.4 Reduced order disK recepte n Cfi 

Eauation (3 10) governs the relation between the actual interactive forces and the actual 
-JTL-*-* "odes a, the interfaces Bt.ua, ions (3 23, - £*£ 

relation between actual degrees of freedom and the reduced order ones. Jo fmd fte reh 
between the reduced order forces f J and the reduced order displacements u ^ 5 5 
the inverse functions for equation (3.23) and (3.35). For equatton . Om • 

overdetetmined by u p, it is possible to find a least square approx, matton of u „ 

error 


Error = 


- d d " d 

U P P.o u 0 


(3 38) 


is minimized. The leas, square approximation of Z l can be mathematically written in a 
simple form 


- d d + - d 

U o = Qp,o u P 


(3 39) 


where Q d p 0 * is the generalized inverse of Q d p , c with the expression 


.1 T „d 


q;; -k, q 


P,0 


Q 


p.o 


(3.40) 


For equatton (3.35), since is.under.de, ermined by f * there are -any posstb 

approximations of f J for every given f t However, by requesting, the nonn of f , «o be 
minimized, i, is possible for us to find a unique approximation off , that retatns the non- - 
equilibrated forces only and neglects the self-equilibrated forces. The express, on 
self-equilibrated approximation is 


d „ d 

f p =\Qp.o 


+i T -d 


f o 


(3.41) 


where Q d p 0 + is defined in equation (3 40) 

Here we are actually proposing an assumption that the self-equilibrated forces a, the 
interfaces are negltgible. This assumption should be valid for the lower system m^e n 
the disk does no, locally deform very much in these cases. It, however, may 
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Combining equations (3 10), (3.39), and (3.41), we can 


have the following equation 


— d 
u O 



(3 42) 



(3 43) 


and 


-d d + 2 6 

<t> j,o = Qp,o ^j-P 


(3.44) 


R d is the reduced order disk receptance matrix associated with the d.sk interfaces and as 
dimension 6-NB by 6-NB * , c is the reduced order interfacia. disp.ac*men« or th^-«h d s 
mode. Equation (3.42) is the reduced order disk receptance equation which g over ns At 
relation between the reduced order displacements u „ and the reduced order fore n 
Member .ha, this equation is the reduced order foim of equation ,3,0). 
equation associated with the disk-blade interfaces. Based on the assumption t he s* 
equilibrated forces at the interfaces are negligible, equation (3.42) approximates 
behavior by 6 degrees of freedom per interface 


3.2 Formulation for the Blades 


3.2.1 Receptance form ulation aiflflls 

To begin with, le, us consider a single blade which is clamped a. its base. The blade base 
moves rigidly so that its motion can , be fully described by 6 degrees of freedom, namely^ 
translations Id 3 torn, ions Le, I ! be the 6-degree-of-freedom 
blade base, where O is a reference point residing at the blade base^ sung ’ 
denote die group of blade nodes which reside a, the blade base and subsenp, a to denote 
dte nodes other than die nodes of group p, the equations of motion for the entire blade 

written as 


C-16 



(3.45) 


M 




-b - b h "* b 
X« = Za + VZ a ,o X 0 


(3 46) 


- b ~b - b 
X p - Qp, 0 X 0 


(3 47) 


where M b C? and K * 5 are the mass, damping, and stiffness matrices associated with 
TZr° nod':: ^ are the displacements of the blade. Q x . *£ 

displacement vector of the group-a nodes assuming .ha, the b.ade follows the base trnd *» 
no. deform a. all ll is the deformation of the group-a nodes relanveto the utf 
configuration f „ e in ' is the external force acting on the blade M „.„Q„.o* • 15 
inertia, force associated with the blade base motion that is introduced when wntmg the 
equations of motion in terms of relative motion The delations of matnces Q^an U „ „ 
are similar to the process described in Section 313; we will not repeat it ero 
bladed disk is operated a. a stationery speed, the steady-state response of the blade 

written as 


-b - b 
X a = U a 


int 


-b - b 
x p = u p 


iai 


“*b - b iQl 

x o = u o e 


(3 48) 
(3 49) 
(3.50) 


Substituting equations (3.48), <3 49), and (3.50) back into equations (3.45), (3.46), and (3.47) 
and eliminating time-dependent terms, we have 


Uu=R b 0 . a iM« 2 <a Mb «.» +, ) Qb -“» 

(3.51) 

- b A b " b 
U p =Qp,o U o 

(3.52) 

b is the receptance matrix associated with the group-a nodes 
L a,a r 

It has the expression 
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(3.53) 




Q 2 M b 


a, a 



-1 


As discussed in Section 3.1.1. calculating the inverse of the dynamic impedance matrix 
l^ > aa -Q 2 M b aa + ifiC b aa ) is undesirable since the matrix tends to be ill-conditioned 
when the excitation frequency Q is close to the natural frequencies of the blade. Similar to the 
way we formulate the receptance equation for the disk, the receptance R a a can be cast in the 
following form by applying the standard modal analysis technique. 






_ . ^ b „b 
+ 2 l Q to j C j 


(3.54) 


where <J> j i0t is the j-th mode for a blade clamped at its base, and m j, £ j, and co j are the modal 
mass, modal damping ratio, and natural frequency of the j-th mode It should be emphasized 
that the degrees of freedom of the group-p nodes do not appear in equation (3.54) since <}) j,p, 
the mode shapes at the blade's base, are zero. 


3.2,2 Force on the blade base 

While equations (3.51) and (3.52) enable us to calculate the motion of the blade for any 
given external excitation f a and blade base motion u G , we have not yet determined the 
resultant force on the blade base The resultant force on the blade base, denoted as f 0 
hereafter, needs to be known in order to assemble the equations for the disk-blade interfaces 
To achieve this, we first consider the blade as an object which can move and deform in space 

** b 

The forces acting on the blade are the excitation force f a and the interactive force at the blade 

-* b) 

interface f G . Since the linear momentum and angular momentum have to be conserved, we can 
write the conservation law, about the reference point O, for small amplitude vibrations, 


I 


- b ^ “* b 
L k = 2, n k e 


int 


(3.55) 


k k 


- b - b 

£ H k +p k* 


- b „ hb — ' b - b\ 

Lk=L' T k + Pk x ^kl 


ittt 


(3.56) 


k k 

b b b 

where L k is the linear momentum^ c^f the k-^h blade node, H k is the angular momentum of the 

k-th node about the node itself, 7t k and x k are the force and the moment acting on the k-th 
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blade node, and p 5 is the position vector from the reference point O to the k-th node By 
defining the following 2 vectors, 


- T 

P = 


-T 

f = 


_ b T b T - b T^ b T - b T b T 

L \ Hi • L k H k • • • L bn H B n 


^b T -b 

n i x 1 


-b T -b T 
71 k t k 


^b T -b T 
71 BN X BN 


(3.57) 


(3.58) 


where BN is the number of blade nodes, equations (3.55) and (3.56) can be rewritten in matrix 
form. 


0 b T - b n bT f b e ifll 

Qo P =Qo f e 


(3.59) 


where Q b 0 is a 6fBN by 6 matrix which can be derived by the same 

SecUon 3.1. Q „ has the function of calculating the resultant momentum about the refere 

point O. 


Now, by definition, the momentum of the blade can be written as 


p b = M b x 


(3.60) 


where M b is the mass matrix of the entire blade. For the steady-state vibration, 


b - b ifll 
x = u e 


(3.61) 


Substituting equations (3.60) and <3 61) into equation (3.59) and eliminating the time- 
dependent terms, we have 


-fi 2 Q b 0 T M b u b = Q b 0 T f b 


(3 62) 


Recall that equations (3.51) and (3.52) describe the displacements of the blade nodes By 
substituting them into equation (3.62) and recognizing that 


f o = Q b p, 0 T f p 


(3 63) 
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we have 




where 


(3 64) 


zfV= - 


0,0 


O' 


Q; M b Q b 0 -O^Q— M D a>a R b a . a M b a>a Q b a ,o 


(3.65) 



,2 


= - ft Qa,o M a, a Ra, a -Qa.o 


(3.66) 


Equation (3.64) gives the expression of f o, the interactive force acting on the blade base, in 
terms of the blade base displacement u 0 and the external excitation f „ The calculation of the 
coefficient matrices^ 0f and H 0A are tedious but straight forward Notice again that the 
blade receptance R a a will be calculated using the blade modes (equation (3.54)) instead of 
directly finding the inverse of the dynamic impedance (equation (3.53)). This remark is 
important since equation (3 54) is efficient and stable when implemented computationally. 


2.2.3 Formulation for all blades 

Equation (3.64) is wntten for the j-th blade Including the subscript “j” it should be 
written as 



(3.66) 


Therefore, for all of the blades, the equation that governs the states of the blade bases is 




where 


(3.67) 



(3.68) 

(3.69) 

(3.70) 


C-20 



< 37,) 

H b 0A = diag[(H b 0 J i ••• Koji ’ ”( H oJnb] (3.72) 

,otice that and ate block-diagonal matrices with diagonal blocks (z?Jj and 


(3.72) 


3.3 Assembling Disk and Blades 


Recall that the purpose of the analysis is to develop a compute " * us 

solve for the steady-state vibration of the blades Equattons q p ). * C ^ motio „ 

calculate the blade vibration provided that the externa. excttanon ^ \ for the blade 

i o are given. Since the external excitation is usually given, our goal 

base motions 

While equation ,3,2) governs the motion of the « * £ " 

interfaces, equation (3.67) governs the motton o te -d-d-b -I equations 

interfaces. Since there are 4 unknown vectors, namely, u o. *. n£-<^ ^ 

We need to introduce two additional equattons to gov raomcntu m 

interfaces. The first of these is Newton’s Interacts Force Law whtch consent 

at the interfaces, 


- d ■* b 

f o = _ f o 


(3.73) 


The second is Hooke’ s Law which provides a constitutive relation a, the interfaces, 

f O = K o,o(u o “ q o) < 3 74) 

l K „ a matrix which represents the stiffness of the disk-blade attachments. Since the 

— ILss is in genera, —-^ 7 ^ 7 ,)^^ 
known quantity. By solving equations (3.42), (3.67), (3.73), l 

can express the 4 unknowns in terms of the external excitation f „ 
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(3.75) 



-b 
fo = 


• + <o(R d o.o + K;'J' 1 H t, 0(I f 0 b 


- d D d 

u o = -R 


0,0 


I + 


^o.o(R d o.o + K- 0 1 J" 1 H ,, oa f 0 b 


- d 

f o = — 


I + Z b „.o(R d o.<,+Ko.J’'H b w f« 


(3 76) 

(3.77) 

(3.78) 


Recall that, equations (3 51) and (3.52) express the blade vibration in terms of the external 
excitation and the blade base motion. Therefore, the vibrations of the blades can be easily 
derived by substituting equation (3.75) into equations (3.51) and (3.52). 


It is important to notice that the most. computationally intensive part in equations (3.75) 


is to calculate the inverse of the matrix 


I + 


Z b o.o(R d o.o+K-J_ 


The dimension of this 


matrix is 6*NB by 6*NB which is significantly smaller than the dimension of a general finite 
element model of a bladed disk For a typical bladed disk with, for instance, 65 blades, the 
finite element model could have thousands of degrees of freedom to be solved for; however, 
the number of the unknown degrees of freedom is only 390 (= 6*65) if we use the reduced 
order model 


Chapter 4: Comparison Of Results 

In Chapter 3, we have developed the algorithm for the reduced order modeling to solve for 
the forced vibration of bladed disks In this chapter, the algorithm will be tested on a 
simplified model problem A code named LMCC (Linear Mistuning Computer Code) has 
been developed based on the algorithm. The results of LMCC will be compared with the 
results of the direct finite element approach to examine the validity of the reduced order 
model 
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4.1 Attachment Stiffness and Damping Ratios 


Before we implement the proposed algorithm, there are still several parameters of the 
reduced order model left to be determined. The first is the attachment stiffness matrix K 0 < 
which simulates the stiffness of the disk-blade attachments (see equation (3.74)). Because of 
the simplified geometry of our test problem, the attachment stiffness is effectively infinite 
Equation (3.74) then implies the continuity in displacements 


- d - b 

U o = u o 


and equations (3.75), (3.76), (3.77), and (3.78) become 


~ b - D d 

u o — - R 


0,0 


1 


-1 b - b 


- b 

f o = 


1 


-1 


H b f b 

H OA * a 


- d D d 

u o = -R 


o.«i +<</<>.<, 


b “* b 

H 0A fa 


- d 

f o = - 


1 + <o« i M 


-1 


b - b 

H 0 A f « 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 


The other parameters of the reduced order model to be determined are the damping ratios 
of the substructures’ modes Since the mechanism of the structural damping is not well- 
realized so far, it is difficult to accurately determine these damping ratios by analytical 
methods However, since the peak amplitude of a slightly damped system is proportional to 
the reciprocal of its damping ratios, we are able to determine these damping ratios by 
numerical experiments. The procedure was used is the following: 

1. Assume the damping ratios of the substructures’ modes are equal to some small value. 
A better choice is the damping ratios of the system modes (if they are available) 
because, based on our experience, the damping ratios of the substructures’ modes are 
in general approximately equal to those of the system modes. 

2. Adjust the damping ratios of the substructures’ modes such that the peak amplitude of 
the tuned response of the reduced order model matches that of a reliable source, for 
example, the physical experiment 
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In the following section, we will proceed to design a test problem to examine the validity 
of the reduced order model 


4.2 Design of the Test Problem 


When designing an appropriate test problem, it is necessary to consider the following 
issues: 6 


1 . 


The test case should represent a realistic object which can vibrate in three-dimensional 
space, pnng-mass models should be avoided because they are in general unrealistic 
Models with beam-type blades will not be considered as appropriate choices either 
since we intend to check the appropriateness of the rigid blade base assumption in 

other words, we want to des.gn a test problem having multiple nodes at each disk- 
blade interface 


2 . 


The test problem should be simple in geometry so that we can 
results based on our physical intuition. 


readily check the 


e whole system should not have too many degrees of freedom so that it can be 
analyzed by the finite element method without costing too much computational time 
or generating the numerical stability problem with the finite element solution 


After considering these requirements, a 
simplified bladed disk with twenty four blades 
(Figure 4.1) was chosen to serve our purpose 
The finite element model is constructed by using 
plate bending elements with three degrees of 
freedom per node, namely, one translation 
perpendicular to the plane and two rotations 
parallel to the plane The thickness of the disk is 
0.8 inch and that of the blades is 0.15 inch. The 
geometry is chosen so that the disk and the 
nominal blade have the natural frequencies that a 
typical turbo-pump stage would have. Both the 
blades and disk were assumed to have material 
properties associated with a super Nickel alloy. 



Figure 4-1 : Geometry of the test 
problem (scales are in inches) 
Thickness of the disk: 0.80 in 
Thickness of the blades. 0. 1 5 in 
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4.3 Tuned System Responses* 


4.3.1 Resonant frequencies and mfldfi StiapgS 

From [7], we know that, if the blades are perfectly tuned, the j-th engine order traveling 
excitation can excite only a certain family of modes. This family of modes has the 
characteristics that the phase difference between the neighboring substructures j <p. 

N is the number of blades and 



(4.6) 


Based on the above knowledge, we can find the resonant frequencies for event family of 
system modes associated with a specific integer j by finding the forced response peaks of the 
bladed disk that result from the j-th engine order excitation 



Figure 4.2 shows the resonant frequencies of the tuned bladed disk predicted by two 
different techniques. The horizontal axis is the integer j mentioned above and the verucal axis 


• The reduced cider model used in this section has infinitely targe anachme « stiffness for disk-blade interfaces 

and experimentally decided damping parameters as discussed in Section 
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is the resonant frequency predicted. The first technique, whose results are deplete y 
diamonds, is the direct finite element approach. Using this method, we first find the modes o 
the entire bladed disk, and, secondly, use the standard Mode Superpos.tion Method to 
calculate the forced responses of the system The second technique, whose results are 
depicted by stars, is the reduced order approach It is apparent that there is reasonably good 
agreement between the frequency predictions of these two approaches. It is also apparent 
that LMCC predicts slightly higher frequencies, probably because the assumption of rigid 
body motion at blades' bases slightly over-constrained the system. 


It is worth noticing that the resonant frequencies of the system follow the blade a one 
frequencies and disk alone frequencies in a particular pattern (Figure 4.3). The blade ataie 
frequencies, which are depicted by the horizontal dashed lines in Figure 4.3 are the " a " ra 
frequency for a blade rigidly constrained at its base. For this particular blade model, the first 
blade alone frequency, 4893 Hz, is associated with the first easewise bending mode he 
second blade alone frequency, 21387 Hz, is associated with the firs, torsional mode, and the 
third blade alone frequency, 30559 Hz, is the second easewise bending mode. There exists 
stiflwise bending modes because the type of finite elements we chose does not have degrees 
freedom in the stiffwise direction The disk alone frequencies, which are depicted by squares, 
are the natural frequencies of the disk with a clamped inner edge and a free outer e ge 



Figure 4-3: Relation between the system frequencies and the 
natural frequencies of the substructures 
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Let us now classify the system modes into four groups as indicated by the four solid lines 
in Figure 4.3 and number these groups from the bottom to the top of the figure. It is observed 
that the group-1 system frequencies follow the natural frequency of the first blade-alone 
bending mode very closely. When the system frequency is close to the blade alone frequency 
and far away from the disk alone frequency, the corresponding system mode has the tendency 
that the disk remains still and the blades vibrate in the associated blade alone mode. Examples 
of this kind of “blade-type” system modes derived from the finite element analysis are shown 
in Figure 4.4 and 4.5 Figure 4.4 shows a system mode with the phase difference between 
blades being 4 Acp. For this particular mode, the blades vibrate predominantly in their first 
bending mode and the disk hardly participates in the vibration though it does have a four- 
nodal-diameter pattern. Figure 4.5 shows a system mode with the same phase difference 
between blades. For this particular mode, the blades vibrate predominantly in their first 
torsional mode and the disk hardly participates 




a) Three-dimensional plot for the mode b) Contour plot for the disk portion 

shape 

Figure 4-4: Mode #9 of the bladed disk (natural frequency = 4682.3 Hz). For this particular 
mode, the disk hardly participates in the vibration. The integer j = 4. 


The group-2 system frequencies follow the disk alone frequencies initially (from j=l to 
j=6) and then follow the natural frequency of the first blade-alone torsional mode (from j=8 to 
j=12). When the system frequency is far away from the blade alone frequency but close to 
the disk alone frequency, the corresponding system mode has the tendency that the disk 
vibrates in the associated disk alone mode and the blades follow the disk without much 
deformation. Figure 4.6 shows an example of this kind of “disk-type” system mode derived 
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from finite element analysis where the phase difference between blades are 4 A<p For this 
particular mode, the blades follow the disk’s vibration passively while the four-nodal-diameter 
disk modes vibrate vigorously The last few system modes in group 2 (from j=8 to j=12) are 
the modes where the disk tends to remain still and the blades tend to vibrate predominantly in 
their first torsional mode. 



Figure 4-5: Mode #57 of the bladed disk (natural frequency = 21260 Hz). For this particular 
mode, the disk hardly participates in the vibration. The integer j = 4. 



Figure 4-6: Mode #33 of the bladed disk (natural frequency = 9714.5 Hz). For this particular 
mode, the disk strongly participates in the vibration. The integer j = 4. 


C-28 






-€- 4-2 0 2 4 € 

a) Three-dimensional plot for the mode b) Contour plot for the disk portion 

shape 

Figure 4-7: Mode #39 of the bladed disk (natural frequency = 19265 Hz) This is a mixed- 

type mode. The integer j = 7. 




c 




a) Three-dimensional plot for the mode b) Contour plot for the disk portion 

shape 

Figure 4-8: Mode #67 of the bladed disk (natural frequency = 27593 Hz). This is a mixed- 

type mode. The integer j = 9. 

In addition to the blade-type and disk-type system modes, there are also the “mixed- 
type” system modes. System modes of this type often have natural frequencies close to or 
far from both blade alone frequencies and disk alone frequencies. They are in transition 
between blade-type and disk-type behavior. Examples of the mixed-type system modes are 
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shown in Figure 4.7, for a group-2 system mode with 7 A<p phase difference between blades, 
and in Figure 4.8, for a group-3 system mode with 9 A<p phase difference between blades. 


4.3.2 Tuned response curves 

After confirming that the reduced order model works reasonably well in predicting the 
system frequencies, we would like to see if the response curves predicted by the reduced 
order model agree well with those predicted by the direct finite element approach For the 
purpose of comparison, we will also present the result predicted by BLDVEB [8], BLDVIB 
is a numerical code for the vibration analysis of bladed disks. It uses a spring-mass model to 
represent the bladed disk 

Figure 4.9 (a), (b), and (c) are the tuned response curves of the blade tips predicted by the 
three different methods In this particular case, the bladed disk is subjected to a fourth engine 
order excitation and the participating system modes are the blade-type modes with blades in 
bending. It appears that the three curves agree well with each other though the response 
predicted by the reduced order model shifts to a slightly higher frequency. Figure 4.10 depicts 
similar plots for a blade-type mode with blades in torsion. It is observed again that responses 
predicted by both LMCC and BLDVIB agree with the response predicted by the direct finite 

element approach 


4.4 Mistimed System Responses* 

In this section, we will compare the mistuned response of the bladed disk predicted by the 
three approaches mentioned in the previous section The bladed disk is mistuned in the sense 
that there are frequency variations from blade to blade. To ensure that the three approaches 
solve for the response of the same mistuned bladed disk, we modified the code of BLDVIB so 
that it can read the mistuned blade alone frequencies from the same source that LMCC does 
For the direct finite element approach, we simulate the blade frequency variations by changing 
the Young’s modulus of each blade. Typical response curves for the blade-type, disk-type, 
and the mixed-type vibrations will be examined 


* The reduced order model used in this section has infinitely large attachment stiffness for disk-blade interfaces 
and experimentally decided damping parameters as discussed in Section 4.1. 
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Figure 9: Tuned response curves using three different 
approaches. System is subjected to a fourth engine order 
excitation. Blades predominantly vibrate in a first bending 

mode. 


Figure 10: Tuned response curves using three different 
approaches System is subjected to a fourth engine order 
excitation. Blades predominantly vibrate in a first torsional 

mode. 
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4.4.1 Blade-tvpe system vibration 

Figure 4.11 (a), (b), and (c) are the typical mistuned response curves of the blade tips 
predicted by the three different approaches. In this case, the bladed disk is subjected to a 
fourth engine order excitation and the participating system modes are basically the blade-type 
modes with blades in bending The tuned response induced by the same excitation was 
discussed in the previous section (Figure 4.9). It appears that the reduced order model and the 
spring-mass model give equally good results 

Figure 4.12 shows the mistuned response curves of the blade tips predicted by the three 
different approaches In this case, the bladed disk is subjected to a fourth engine order 
excitation and the participating system modes are the blade-type modes with blades in 
torsion. The response of the tuned stage was discussed in the previous section (Figure 4.10). 
We can see that LMCC works much better than BLDVIB, though they predict equally well 
the results for the tuned stage One of the possible explanation for the failure of BLDVIB in 
this particular case is that the spring-mass model does not account for the effects of three- 
dimensional motion and multiple node attachment at the disk-blade interfaces. By observing 
Figure 4. 12(c), we also find that the blades do not interact as much as they should This leads 
to another possible reason for the failure of the spring-mass model that they can not account 
for the complicated blade-to-blade coupling effect while the blades vibrate predominantly in 
their first torsional mode 


4 A. 2 Disk-type system vibration 

Figure 4.13 shows the comparison of the mistuned results predicted by the direct finite 
element approach and the reduced order model. For this particular case, the bladed disk is 
subjected to a fourth engine order excitation and the disk participates vigorously in the 
vibration since the natural frequency of the four-nodal-diameter disk modes is very close to 
the system frequency. It appears the result predicted by the reduced order model agrees well 
with the result predicted by the direct finite element approach. It is interesting to observe 
that the blade mistuning does not affect the disk-type vibration very much. 

4.4.3 Mixed-type system vibration 

Figure 4.14 and Figure 4.15 are the mistuned responses for the mixed-type system 
vibration The vibration is mixed in the sense that both the disk and the blades participate in 
the vibration and it is hard to tell which blade alone mode the blades vibrate in. Figure 4. 14 
shows the response curves for the mistuned bladed disk under a seventh engine order 
excitation. Figure 4.1 5 shows the response curves for a mistuned bladed disk under a ninth 
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Figure 1 1 : Mistimed response curves predicted 
by three different approaches System is subjected 
to a fourth engine order excitation. Blades 
predominantly vibrate in a first bending mode. 


Figure 12: Mistimed response curves predicted 
by three different approaches. System is subjected 
to a fourth engine order excitation. Blades 
predominantly vibrate in a first torsional mode. 
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Amplitudes of BlnHe Tips f!ni 


a) Direct finite element b) Reduced order model 

Figure 4-13: Mistuned response curves predicted by different approaches. System is subjected to a 
fourth engine order excitation The disk strongly participates in the vibration 



a) Direct finite element b) Reduced order model 

Figure 4-14: Mistuned response curves predicted by different approaches. System is subjected to a 
seventh engine order excitation Both the disk and blades participate in the vibration 




Figure 4-15: Mistuned response curves predicted by different approaches. System is subjected to a 
ninth engine order excitation. Both the disk and blades participate in the vibration. 
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engine order excitation. It is apparent that the curves predicted by the reduced order model 
agree very well with the curves predicted by the direct finite element approach. 

In the last two sections, we observed that the reduced order model can accurately predict 
the response of both tuned and mistuned bladed disks. The spring-mass model used by 
BLDVIB failed to model the coupling effect between blades while the blades vibrate 
predominantly in their first torsional mode In the next section we will compare the execution 
time taken by the direct finite element approach and the reduced order approach 


4.5 Execution Time Comparison 


Figure 4.16 shows the execution time comparison between the direct finite element 
approach and the reduced order approach when using a supercomputer. 



Figure 4 16: Execution time comparison between the direct finite element 
approach and the reduced order approach. 


The data are calculated based on the estimation for a typical bladed disk as following: 


Number of blades 50 

Number of degrees of freedom of a blade 1 5,000 

Number of degrees of freedom of a disk sector 10,000 

Number of degrees of freedom of a bladed disk 1,250,000 
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The estimated CPU time for solving for 100 system modes is 1,000 CPU minutes per bladed 
disk. The time for calculating the forced response using Mode Superposition Methods is 
approximate 0.02 minute for every 100 excitation frequencies. Therefore, we can write down 
the following equation for the relation between the execution time needed by using the direct 
finite element approach and the number of bladed disks analyzed, 

T fem = 1000.02 N (4.7) 

where N is the number of bladed disks. 

If we use the reduced order approach for the blade mistuning problem, the CPU time for 
solving for the modes of a nominal blade is approximate 15 minutes and that for the disk 
modes is approximate 250 minutes provided that the cyclic symmetric boundary conditions 
are used to improve the computational efficiency. The CPU time needed to process the 
output of the finite element analysis for the use of LMCC is approximate 15 minutes And, 
for every 100 excitation frequencies, LMCC will take approximate 5 minutes to calculate the 
forced response Therefore, we can write down the equation for the relation between the 
execution time needed by using the reduced order approach and the number of bladed disks 
analyzed. 


jlmcc _ 2g 0 + 5 N ( 4 g) 

An immediate comparison we can obtain from observing equations (4.7) and (4.8) is that, for 
every additional bladed disk, the direct finite element approach needs 1000 more minutes but 
LMCC (the reduced order approach) needs only 5 more minutes. Figure 4.16 shows that, for 
a typical mistuning analysis with 100 bladed disks analyzed, the direct finite element 
approach needs 100,002 minutes and the reduced order approach needs only 780 minutes. 


Chapter 5: Conclusions And Suggestions 


■ 5.1 Conclusions ■ 

Mistuned vibration of bladed disks has drawn much attention since it can result in large 
blade-to-blade amplitude variations and the high response blades can fail from high cycle 
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fatigue. Due to the severe operating conditions that modem aircraft are subjected to, accurate 
predictions of the vibration of mistuned bladed disks become necessary. Using finite elements 
to model a bladed disk is not a feasible approach for a thorough mistuning analysis since the 
geometry of a bladed disk is so complex that a structurally accurate finite element model 
would have too many degrees of freedom Using spring-mass models to represent a bladed 
disk results in low degree-of-freedom formulations Though they can be efficiently and 
reliably solved by computers, the ad hoc way in which they choose system parameters and 
their questionable accuracy are undesirable 

To develop a structurally accurate and computationally efficient model for the bladed 
disks, we first considered the receptance method Though the receptance method can 
significantly reduce the number of degrees of freedom of the system, the resultant equations 
can still have too many degrees of freedom since the disk and the blade usually contact at 
multiple points Its necessity of using free-free blade modes is especially undesirable since, in 
general, blades tend to vibrate like clamped-free modes and much of the work done on the 
blades utilizes cantilever blades. 

The reduced order modeling follows the idea of the receptance method but further 
assumes that the blade bases undergo rigid body motion. The blade vibration is considered as 
a combination of blade base motion and cantilever blade modes. By finding the six equivalent 
degrees of freedom for the disk vibration at each disk-blade interface, we are able to reduce the 
number of degrees of freedom of the whole structure down to six times the number of blades 
With this significant reduction in the number of unknowns, the reduced order model results in 
a formulation that can be efficiently solved by computers. Along with its low computational 
cost, the reduced order model is structurally accurate and provides a straight forward 
approach for modeling the system since it directly takes the results from the finite element 
analysis of a tuned system without ad hoc assumptions. 

The validity of the reduced order model was examined by applying its algorithm to solve 
for the forced response of a carefully designed test problem. The test problem is designed so 
that it represents a realistic object that vibrates in three-dimensional space and has multiple 
nodes at each disk-blade interface For the purpose of comparing the reduced order model 
with other reliable techniques, the geometry of the test problem is relatively simple so that it 
can be efficiently solved by the direct finite element approach. It is shown in Chapter 4 that 
the reduced order model can predict the tuned system frequencies reasonably well. By 
comparing the results from the different approaches, it appears that the reduced order model 
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is able to capture the characteristics of the mistuned bladed disk. It predicted response curves 
that are almost identical to the results of the direct finite element approach for the more than 
one hundred cases examined, of which the cases presented in Chapter 4 were representative. 
The results predicted by a typical spring-mass model are also examined. It appears that the 
spring-mass model predicts mistuned responses well when the blades vibrate in bending 
modes but predicts poorly when the blades vibrate in high frequency torsional modes. 

Based on the investigations we have done, we can conclude that the reduced order model 
is a valuable model for the vibration analysis of mistuned bladed disks It is valuable for the 
following reasons: 

1. It is structurally accurate and can be used for bladed disks of various kinds of 
geometry because it utilizes the output of general finite element models instead of 
being restricted to springs and masses or plates and beams. 

2. It is straight forward in the sense of modeling since the system parameters are directly 
taken from the output of the finite element analysis instead of making ad hoc 
assumptions 

3. It is computationally efficient because the number of degrees of freedom of its 
formulation is only six times the number of blades. This number is significantly 
smaller than that of a finite element model and of the same order of that of a more 
complicated spring-mass model, for example, [9], 


5.2 Suggestions 

One of the subjects of bladed disk vibration under intensive study nowadays is the friction 
related mistuning problem. Blade-to-blade and blade-to-ground friction dampers are integrated 
into the bladed disk assemblies to dissipate energy so that the peak amplitudes of blade 
vibration can be reduced. Since the nonlinear friction forces are sensitive to the displacements 
and, sometimes, to the velocities of the blades, it is necessary for the analysts to have a 
structurally accurate model. Much of the work has been done on this subject is based on 
spring-mass models. With proper tuning, spring-mass models often predict reasonably 
accurate responses for a linear system, however, their ability in predicting the response of a 
nonlinear system is questionable since friction forces can considerably change the mode 
shapes of the blades and this change can hardly be modeled by simple spring-mass systems. 
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It is, therefore, reasonable to use the reduced order model for the analysis of friction 
related mistuning problems. Because of the computational efficiency of the reduced order 
model, we can easily derive the receptance for the friction interfaces of the bladed disk It is 
then possible to use this structurally accurate receptance to solve for the friction problem 
either by time integration methods or by approximation techniques. 
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